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Abstract. Photon coincidence spectroscopy is a promising technique for probing the 
nonlinear regime of cavity quantum electrodynamics in the optical domain, however 
its accuracy is mitigated by two factors: higher-order photon correlations, which 
contribute to an enhanced pair count rate, and non-simultaneity of emitted photon 
pairs from the optical cavity. We show that the technique of photon coincidence 
spectroscopy is effective in the presence of these effects if the quantitative predictions 
are adjusted to include non-simultaneity and higher-order correlations. 
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1. Introduction 

Photon coincidence spectroscopy (PCS) provides an accurate means for detecting 
quantum field effects in cavity quantum electrodynamics (CQED) (Carmichael 
et al. f996, Sanders et al. f997, Horvath et al. f999, Horvath and Sanders 2001). In 
applications of CQED to quantum information, where entanglement between the cavity 
field and the internal degrees of freedom of the atom passing through the cavity is 
central to the scheme (Turchette et al. 1995), PCS offers an accurate method for probing 
the quantum field in the cavity. In two-photon coincidence spectroscopy (2PCS), the 
atom-cavity system is driven by a bichromatic field (as for nonlinear spectroscopy; 
see Thompson et al. 1998), and a detector measures the two-photon coincidence 
rate (2PCR) in the cavity output field. Two-photon spectral peaks are expected for 
particular sum frequencies of the bichromatic driving field (one frequency is fixed and 
the other scanned across some range of frequencies). Certain two-photon spectral peaks 
would occur only if the quantum description of the intracavity field is valid and not occur 
if the intracavity field can be described by a semiclassical field theory. 

In order to test the quantum description of the intracavity field, accuracy in 
predicting the locations of the two-photon spectral peaks is very important. The height 
or the width of the two-photon spectral peaks does not need to be determined by 
theory accurately, but the centre of the peak, as a function of scanning field frequency, 
is important. We consider here two effects on the shifts of these peaks: higher-order 
photon correlations that contribute to the 2PCR and the non-simultaneity of the photon 
pairs as they exit the cavity. These two cases contribute to small shifts of the 2PCR 
peaks but must be accounted for in establishing accurate predictions by 2PCS. 

2. Master equation and two— photon count rates 

In the electric-dipole and rotating-wave approximations, the Jaynes-Cummings (JC) 
Hamiltonian for the two-level atom (2LA) coupled to a single mode is (Jaynes and 
Cummings 1963) 

H(g) = lo(g z + a) a) + ig(a)a_ — aa + ), (1) 

with g a position-dependent dipole coupling strength, a and the annihilation and 
creation operators for photons in the cavity field, o~ + , o~_, and a z the 2LA raising, 
lowering and inversion operators, respectively, and H — 1. Single-atom cavity QED 
in the optical domain (Hood et al. 2000, Pinkse et al. 2000) validates the single- 
atom approximation of Eq. ([!]). Provided that the atoms move sufficiently slowly 
through the cavity (Carmichael et al. 1996, Sanders et al. 1997), and the position is 
randomly varying, the coupling strength g is also random: hence, a coupling strength 
distribution P(g) can be constructed (Sanders et al. 1997), and we assume the P(g) 
depicted in Fig. 5 of Sanders et al. 1997, which applies for a single TEM o cavity mode 
and atoms travelling transverse to the cavity axis. The atoms first pass through a 
rectangular mask, centred at an antinode with size u>o x A/10 (with wq the cavity mode 
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waist and A the optical wavelength). The atoms then traverse the cavity after passing 
through the mask. We restrict Fg max < g < g max for g max the coupling strength at an 
antinode along the cavity longitudinal axis and F an effective cut-off term. 

The JC spectrum for the Hamiltonian ([I]) is depicted in Fig. |I] and the 'dressed 
states' of the atom-cavity system are designated as the lowest-energy state |0) = 
|0)cav ® | g) atom = 1 0, g) , and, for n a positive integer, the 'excited' state couplets \n)± = 
i/y/2(\n — 1, e) ± i \n, g)), with \n) the Fock state of the cavity mode and |g) (|e)) the 
ground (excited) state of the 2LA. Averaging over P(g) yields inhomogeneous spectral 
broadening (due to atomic position variability). 

Two-photon excitation is provided by driving the atom directly, as it traverses 
the cavity, with a bichromatic field S(t) = £ie~ tuJlt + E 2 e~ luJ2t . The driving-field 
frequency uj\ is fixed and resonantly excites |0) < — > for the subensemble g = 
gf = uj — uji corresponding to P(g) = 5(g — (?/). The scanning-field frequency u 2 
excites |1)_ < — > \2)± for P(g) = 8(g — gf). Enhanced rates of photon pair detection are 
then sought as the scanning frequency u 2 is varied such that u)\ + uj 2 is resonant with 
some transition |0) « — > \2)± as depicted in Fig. [1| 

The master equation for this system (Sanders et al. 1997) can be expressed as 
p = Cp for £ = C cS + V + J, i.e. a sum of a Liouville operator C eS , an explicit time- 
dependent Liouville operator T> and a 'jump' term J . We introduce 5 = uo 2 — u)\ and 
work in the rotating picture with respect to the driving-field component 

For = ig(a)a^ — aa + ), the effective Hamiltonian is 

H eS (g, Ei) = (u - ui) (a z + a) a) + E(g) + T(Ei) 

— iKaJa — i(j/2)a+cr-, (2) 

with 2k the decay rate for the cavity, 7 the inhibited spontaneous emission rate for 
the partially confined atom passing through the cavity, T(£i) = i£i(a + — a J) a 
monochromatic 2LA driving term, and 

£es($,€i)p=-i H eS (g,Ei)p- pHl B (g,Ei) . (3) 

The jump term superoperator J is given by J p = 2napa' f + ■ya-pa + , and Vp = 
-i [r(£ 2 e- iSt ),p\. 

Solving the master equation for p, and accounting for atomic position variability, 
yields a solution 

p= / P(g)p(g)dg. (4) 

For a bichromatic driving field, the density matrix does not settle to a steady state 

value, but, in the long-time limit t — > 00, the Bloch function expansion is 

00 

lim p(t) = Y j p m e im5t , (5) 

m=0 

with p m time-independent matrices. As the photocount integration time is expected to 
be long compared to the detuning 5, it is reasonable to assume that rapidly oscillating 
terms average out and therefore approximate p(g) by truncating the expansion (|5|). 
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An experimental signature for entanglement is the 2PCR, which is given by the 
expectation value 

w (2) M,£i) = (2k) 2 (: n(t )n(t + t) :) , (6) 

for n(t) = a^(t)a(t) and ': :' referring to normal ordering. The coefficient (2k) 2 accounts 
for the photon pair rate of emission from the cavity. The expectation value corresponds 
to a trace over p in Eq. (^) for some choice of P(g)- The 2PCR is of interest in the 
asymptotically large time t limit, where transients disappear and t is unimportant: 
hence the t dependence on the left-hand side of Eq. (||) is ignored. The 2PCR is 
obtained by collecting photons over some duration, which we call the window time r w . 
Thus, the actual 2PCR is an integral, or a double integral, over t of uA 2 )(i, 5, g), yielding 
an effective two-photon coincidence rate of w^(5, g) for a window time of r w . The 
details of the window time, and the choice of performing a single or a double integral, 
are fully discussed in Sec. |J For this Section, we assume that the window time is 
sufficiently small that (Sanders et al. 1997) 

w^(S,S 1 ) = (2K) 2 (:n 2 :). (7) 

For 2PCS, the major obstacle to detect entanglement is the background 2PCR due 
to two photons of angular frequency u 2 contributing to excitation to the second couplet, 
generally by off-resonant, but nonnegligible, transitions. The method for overcoming 
this problem is called 'background subtraction'. In this technique, the experiment is 
performed twice, once with a bichromatic driving field to obtain w^ 2 \5, £\) and again 
but with a monochromatic driving field at frequency u>2 (i.e. S± = 0). The difference 
between the two 2PCRs, for the cases that S\ ^ and Si = 0, is given by 

A^,^) = w^(5,S 1 )-w^(5,S 1 = 0), (8) 

which we refer to as the 'difference-2PCR'. The 2PCR with and without background 
subtraction is depicted in Fig. |, and the importance of background subtraction is 
apparent for the two-photon spectral peaks in the domain \5\ < 1 for 5 = (u>2 —uj)/(uj — 
wi) and P(g) = 5(g - g f ). 

3. Multi— photon contributions 

The emphasis thus far has been on 2PCS, with its two-photon decay process leading 
to an enhanced 2PCR. Although two-photon events dominate the dynamics, multi- 
photon effects cannot be ignored. In fact, multiphoton contributions can be used to 
extract signatures of transitions from higher-order couplets in the JC ladder. This 
method is known as multi-photon coincidence spectroscopy (Horvath et al. 1999). 

In principle multi-photon emissions from the cavity can be distinguished from two- 
photon emissions, but, in practice, a photodetection system located outside the cavity, 
which is designed to detect two coincident photons passing through the cavity mirror, 
cannot distinguish between two coincident photons and the rarer case of N coincident 
photons with N > 2. In addition to the case of multi-photon emissions in the cavity 
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output field, there are also cases where multiple photons simultaneously exit the cavity 
mirror and exit the side of the cavity. Hence, the two-photon coincidence rate includes 
contributions from multiple photons emitted from the cavity as well as two photons 
emitted by the cavity and at least one photon leaving the side of the cavity. 



ir™(t, 5,S l ) = J2 a k( :n - n (*l) ■ II Q k s -M : ) , (10) 



3.1. Correlation Functions 

In order to treat multiphoton emissions, the correlation function w^ 2 \t, 5, £\) in Eq. (|^), 
must be generalized to the multiphoton case w^(t,S,£i), for t > 2 the number of 
photons emitted. These photoemissions take place at random times 

t={t u t 2 ,...,t l „ 1 ), (9) 

with t a vector consisting of the photoemission times. Only I — 1 photoemission times 
are required for I photoemissions, as each photoemission time is determined with respect 
to the time of the first photoemission which can be fixed at to = 0. Thus, the higher 
contributions to the 2PCR are given by Y1T=2 w ^{t, 5, Si). 

A compact notation for is possible if we introduce a single notation for the 
number operator n corresponding to the cavity field and for the atomic inversion 
operator cr + cr_. We introduce the operator 0j, with i £ {0, 1}, and 6 = n, Oi = a + a_. 
By introducing the length I — 2 vector k, such that each hi £ {0, 1} Vi, we obtain the 
compact expression 

/ £ ~ 1 

k \ s = 2 

with {afc} a set of real constant coefficients. The sum is over all 2 distinct vectors 

ke {(o,o,. ..,o,o,o) ,(o,o,. ..,o,o,i), 

(0,0,... ,0,1,0), ...,(1,1,... ,1,1,1)}. (11) 

The operators n and a + a_ correspond to quantum numbers inside the cavity: the 
corresponding extra-cavity terms are 2nn and 7a + a_. These coefficients correspond to 
the rates for the quantum to leave through the cavity mirror and through the side of 
the cavity, respectively. For each n, we introduce the multiplicative factor of 2k, and 
the coefficient 7 applies for each <7 + cx_ in the product. Thus, 

e-2 

a k = (2^y- K 1 K ioiK = J2 k i- ( 12 ) 

i=i 

For example, let us consider £ = 4. In this case 

w (4) ((M 2 ,*3),<y,£i) = (2K) 2 (:n-n(t 1 ) ■ [(2K) 2 e (t 2 )e (t3) 

+ 2^Q {t 2 )Q l {t 3 ) + 2K 1 Q 1 {t 2 )Q {t 3 ) 
+ 7 2 e 1 (t 2 )6 1 (t 3 )] :)• (13) 

We observe that is the correlation function for two photons at the initial time 
to and a later time ti, plus the set of all photons emitted from the cavity or out the side 
at the other times {t 2 , t%, . . . , U-i}. 
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In Sec. f| we consider the window time r w in detail, but here we make the small 
window time assumption kt w <C 1. Hence the t dependence is negligible, and 

n-2 

w (n ^ 2) (^,5) = ^(2K) n - m 7 m (a™a tn - m a n - m a™). (14) 

m=0 

Expression QI^D is valid provided that the window time is small, as we assume in this 
Section. 

3.2. Three-photon coincidence rate 

In Fig. ^| the three-photon coincidence rate (3PCR) is shown for the cases (a^ 3 a 3 ) 
and (V_at 2 a 2 cr + ) with fixed coupling strength g = Qf = 9k. Vertical lines depict the 
seven values of 5 which would lead to on- or off-resonant excitation pathways which, 
after the third photon absorption, are on-resonance with either |3)_ or |3)+. These 
vertical lines serve as a guide to where peaks in the 3PCR may be expected. The values 
of S for which these vertical lines occur are given by the resonance conditions 

3u — V3g = 3uii, (15) 

3u ± V3g = 2uj\ + U2, 8 = 2 ± V3g/gf, (16) 

3u±V3g = u 1 + 2u 2 , 8 = (g f ± y/Sg) /2g f , (17) 

3cu±V3g = 3u 2 , 8 = ±g/VSg f . (18) 

The 3PCR peaks are observed near these vertical lines, although small shifts of the peaks 
are evident. The reason for these peak shifts is the influence of competition between 
excitation pathway, as discussed in Horvath and Sanders 2001 for 2PCS. These shifts 
in peak positions are small and do not affect the analysis here. The total 3PCR is the 
weighted sum of the full and dotted curves of Fig. |3|: (2k) 2 (2k (a* 3 a 3 ) +7 (cr + a^ 2 a 2 a_)). 

3.3. Analysis of the three-photon count rate 

In order to understand the detailed structure of Fig. |3|, we modify the Liouvillean 
superoperator £ by artificially eliminating particular driving terms responsible for 
certain peaks. 

Although a thorough examination of the peak structure in Fig. |3| requires an analysis 
of both rates (a^ 3 a 3 ) and (u+a^aV-), an analysis of one of the two graphs is sufficient. 
The rationale for the peak structure is readily extended to explain the similar peak 
structure observed for the other rate. Without loss of generality, we choose, in this 
Section, to study in detail the peak structure for the 3PCR (a^ 3 a 3 ). This analysis 
of the 3PCR structure is important in Subsection |3.4], as the 3PCR is responsible for 



perturbative corrections to the desired 2PCR, which provides the experimental signature 
of quantum field effects. In addition, the 3PCR structure could be used to identify 
further signatures of quantum field effects compared to the 2PCR signature alone. 

The isolation of specific transitions is obtained by the following procedure. For 
example, let us consider the influence on the 3PCR of the |0) < — ► |1)_ transition. We 
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can write the effective Hamiltonian (0) as a matrix in the dressed-state basis. To isolate 
this influence, we can set the matrix elements (0| Y(£i) |1)_ and T(£i) |0) to be 
zero where T(£i) is the driving term of Eq. (0). In addition the matrix elements of the 
driving term (0| T(£ 2 exp(— iSt)) |1)_ and its complex conjugate can both be set to zero. 
The jump term J is not modified because only the driving terms and their effects are 
of concern in this analysis. 

The result of setting the matrix element (0| T(£i) |1)_ and its complex conjugate to 
zero for calculating the rate ^a^a 3 ) is depicted in Fig. f|(a). Similarly Fig. f|(b) shows 
the graph for (0| T(£ 2 exp(— i5t)) |1)_ and its conjugate set to zero. We observe that 
eliminating these two matrix elements causes a dramatic reduction of neighbouring peaks 
at both 5 = — (v^2 — l), which is depicted as peak ii, and 5 = — (v3— l) /2, which 
is depicted as peak Hi in Fig. 0. The reason for this reduction is that U\ corresponds 
to the |0) < — > |1)_ transition. To understand the next step we refer to Fig. 0(b) 

which depicts ^a^a 3 ) with £ 2 = for the |1)_ < > |2)_ transition. In this figure it is 

clear that the photon of frequency uu 2 dominates the |1)_ < — > |2)_ transition. For the 
peak 5 = —(\^2 — 1), this transition is on resonance, but, for the peak 5 = — (y/3 — l)/2, 
the transition is slightly off-resonant. 

The final contribution to the rate appears in Fig. 0(b) where the contribution of 
an uo 2 photon to the peaks is apparent. However, the transition |2)_ < — > |3)_ for 
the peak 5 = —(V2 — 1) is not resonant, but the transition is resonant for the peak 
at 5 = — (y/3 — l)/2. From the three figures [§(a), 0(b) and 0(b), we observe that the 
peak rate at 5 — — (v^2 — 1) and 5 = — (\/3~ — l)/2 is overwhelmingly due to a sequence 
of an uji photon followed by two u 2 photon absorptions. By removing one of these three 
photons artificially, via setting the corresponding matrix elements of the Liouvillean to 
zero, the peaks at 5 — —(\^2 — 1) and 5 = — (v^— l)/2 are dramatically decreased. The 
excitation pathway for the peak at 5 — — (v3 — l)/2 may be understood in a simple 
way as follows. Once again setting E\ = for the |0) < — > transition reduces the 
peak, not only at 5 = — (y/2 — 1), but also for its close neighbour 5 = —(a/3 — l)/2. 

Let us now consider the increase in peak v of Fig. f|(a). This increase is due to a 
decrease in competition which arises by eliminating the |0) < — ► |1)_ transition due to 
an ui photon. This competition is clear from the following analysis. From Fig. |](d), we 
see that peak v is reduced by eliminating the uo 2 photon which drives the |0) < — > |1) + 
transition. Then Fig. 0(h) shows that peak v is also reduced by eliminating the u> 2 
photon-induced |1) + « — > |2) + transition and, finally, the uj 2 photon is also responsible 

for the |2) + < ► |3) + transition as shown in Fig. 0(h). Therefore, removing the iO\ 

photon-induced |0) < >■ |1)_ transition allows a greater population to occur in the |1) + 

level in Fig. ^(d) and then, via two more to^induced transition, a greater rate (a^ 3 a 3 ). 

In Fig. [|(a) we observe that peak vi is completely eliminated by setting E\ = for 
the |0) « — > transition. We can see from Figs. 0(d) and 0(h) that an uji « — > uj 2 < — ► 
uj 2 pathway exists for |0) < — ► |1)_ < — > |2) + < — > |3) + , which is responsible for peak vi 
in the rate (a^ 3 a 3 ). A second, less important, pathway is evident in Figs. 0(b) and 0(d) 
as the ui < — > u 2 * — > u 2 induced pathway |0) < — > |1)_ < — > |2)_ < — >■ |3) + . Thus, 
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there are two contributing pathways, one via |2)_ and the other via |2) + , which are both 
prevented by eliminating the u>i photon-induced |0) « — ► |1)_ transition. Therefore, the 
peak in Fig. |(a) is reduced for the condition that the |0) < — > is suppressed. 

We observe in Fig. [3] that peak iv is suppressed for peak 5 = 2 — y/3. This involves 
a resonant excitation to |1)_ by an uj\ photon followed by an 0J2 photon exciting to |2)_ 
with a detuning of (l+y/2— Vtygf = 0.68<?/ followed by an uj\ resonant excitation to |3)_. 
Another dominant pathway involves a resonant excitation to |1)_ by an uj\ photon 
followed by another u)\ photon excitation to |2)_ with detuning (2 — y/2)gj = 0.59g/, 
followed by a resonant excitation to |3)_. 

We do not observe a peak at 5 = 2 — y/3 because the off-resonant excitations are 
small in comparison to excitations near the two peaks Hi and v. In contrast, peak viii 
also involves a nonresonant pathway which is similar. This pathway involves a resonant 
excitation to |1)_ by an u\ photon followed by another uj\ photon excitation to |2)_ 
with detuning (2 — y/2)gf = 0.59gj, followed by a resonant excitation to |3) + . Peak viii 
is visible because there is no major peak nearby, but its height is less than the valley 
height at iv. 

To complete the analysis of 3PCR, we consider the intriguing bump on the 
graph in Fig. |3] near 5 = —0.7. There is no vertical line in Fig. |3| at this bump 
because this value of S does not correspond to any resonant excitation pathway. From 
Figs. f|(b), ||(b) and ^|(b) we observe that there is a 3u>2 induced excitation pathway 
along |0) < — ► |1)_ < — > |2)_ « — ► |3)_. The frequency for the |0) < — ► |3)_ transition 
is 3co>2 = 3to — y/3gj which corresponds to 0J2 = co — gj/ y/3. This value of U2 is equivalent 
to setting 5 = —1/ y/3, which corresponds to the solid vertical line i in Fig. ^ The reason 
for the bump occurring at 5 = —0.7 rather than at 5 — — l/y/3 is shown in Fig. [?| and 
explained below. 

In Fig. [7] each excitation pathway involves nonresonant transitions. In each of 
the three cases, the transition |0) « — > |1)_ occurs on resonance (5 = —1) or slightly 
off resonance (5 = —l/y/2 and 5 = —l/y/3). The 5 = —1 case does not generate a 
large population of the |3)_ level, because of the significantly off resonant excitations 
to |2)_ and to |3)_. Similarly, the 5 = — l/y/3 case also does not significantly 
populate |3)_ because of the large detuning from |2)_ compared to the homogeneous 
linewidth 3k + 7/2. The dominant pathway is thus 5 = —l/y/2 = —0.7, for which the 
excitation is close to resonance for |1)_ and |3)_ and exactly on resonance for |2)_: 
hence the observed bump at 8 = —0.7. 

In Figs. ||(c) and (d) we observe that the bump at 5 = —0.7 is reduced by 
competition with the |2)_ < — > |3) + pathway. By setting either S\ = or £ 2 = 
for this transition, the bump is increased. Thus, excitation along the |0) < — > \1)_ < — > 
\2)_ < — > |3) + pathway diminishes the 3PCR, and this decrease is more pronounced 
for 3to> 2 > 3c<j — y/3gf. Therefore, a consequence of this competition between excitation 
pathways to |3)_ and |3)+ is an apparent shift of the peak responsible for the bump to 
lower 5. 

Figures provide sufficient information to extend the analysis above and identify 
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specific excitation pathways to determine which contributions dominate the 3PCR graph 
of Fig. |3|. The cases discussed above are illustrative examples. Also a similar analysis 
can be applied to understand the structure of the 3PCR {o- + a) 2 a 2 aJ) , but we do not 
perform this analysis here. 

3.4- Difference two-photon count rate 

The background subtraction discussed in Section |2] is relevant to multiphoton 
contributions. The difference-2PCR is constructed in the same way as before, except 
that the multi-PCR is taken into account: 

A (n) (5,£i) =w i - n) {5,£ 1 ) -w (n) (S,£ 1 = 0). (19) 

The 2PCR contains contributions from all higher-order terms; thus, when one of the two 
chromatic components in the bichromatic driving field is shut off, the difference-PCR 
includes subtraction of all higher-order correlation functions. 

It is convenient and reasonable to neglect higher-order multiphoton contributions 
to the 2PCR and retain only the term of interest and the lowest-order correction 
term w^ 3 \ Thus, the 2PCR in the perturbative limit can be approximated by +w^ 3 \ 

The difference-2PCR is depicted in Fig. |8] as a plot of A^ 2 ^ versus 5, including the 
3PCR contribution. After background subtraction, we observe 2PCR peaks for 5 = 
1 — y/2, 5 = y/2 — 1 and 5 = v2 + l, depicted as peaks II, III, and V, respectively. These 
three peaks have been identified and discussed in detail in Carmichael et al. 1996 and 
in Sanders et al. 1997. The important point in Fig. [5] is that the 2PCR without 3PCR 
corrections, in Sanders et al. 1997, corresponds closely to the corrected version here. 
That is, the solid and dotted lines in Fig. [S] are quite similar. Therefore, 3PCR would 
not destroy the desired peak in 2PCR. However, three-photon effects are manifested on 
the 2PCR peak in Fig. || and this extra structure is of interest. This noticeable extra 
structure occurs for values of 5 corresponding to the vertical lines I, IV and VI. 

The extra structure due to 3PCR is not surprising. This is because 2PCR requires 
excitation to \2)± with zero detuning in order to observe the best possible 2PCR. This 
condition of excitation to \2)± is also the prerequisite for observing 3PCR, as resonant 
excitation to |2)± ensures a nonnegligible population of |3)±. 

A particularly significant 3PCR effect is the enhancement of the peak at 5 = 
—(y/2 — 1) and a slight shift in the peak position. In Fig. [8] we can compute the 2PCR 
with and without the 3PCR correction. The 3PCR correction term is significant. This 
3PCR contribution applies for values of S which contribute to peaks ii and Hi in Fig. [| 
Although Fig. [3] applies only for g = gf, this represents the dominant contribution 
after background subtraction. Therefore, Fig. |^ is a useful guide for understanding the 
enhancement of peak II in Fig. 

Minor peaks occur at 5 values for vertical lines IV and VI in Fig. ^. The origin 
of the peak IV is the same pathway responsible for peak vi in Fig. [| Similarly, peak 
VI corresponds to peak viii in Fig. [|. The peaks in Fig. |] enable us to understand the 
extra structure of the 2PCR in Fig. ||. A quantum trajectory simulation of the 2PCR, 
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which included the 3PCR corrections, multiple atom effects and finite window times, 
is shown in Fig. 5 of Carmichael et al. 1996. Peak IV is slightly discernible in Fig. 5 
of Carmichael et al. 1996 although subject to a small shift in the peak position (shifts 
in peak position have been discussed in Horvath and Sanders 2001. The information 
in Subsection [O] enables a complete understanding of the structure in Fig. [8] and, in 
practice, the enhanced 2PCR in Fig. 5 of Carmichael et al. 1996. 



4. Non— simultaneity of photodetections 

In the following analysis we investigate the appropriate window time for counting photon 
pairs. The higher-order correlations can be ignored in determining the optimal window 
time r w as the correction for this optimal time will be small. Photon pair emissions, 
which arise due to de-excitation from the |2)± level to the ground state of the JC ladder, 
are not emitted from the cavity simultaneously due to the randomness of photoemission 
due to the cavity linewidth. The detection of a photon pair thus depends on identifying 
a window time r w such that, for two photons detected with temporal separation t < r w , 
the two photons are deemed to be members of a pair, lit > r w , the photons are deemed 
to be independent single photons (not members of a pair). Here we determine the 
optimal window time r opt which maximises the counting rate of genuine pairs relative 
to the rate of false pair counts. 



4-1. The two-photon count rate 

The 2PCR can be obtained in more than one way. Ideally one would have a 
perfectly efficient photodetector that detects all photons leaving one cavity mirror. 
The photodetector would then provide a complete record of photon emissions from 
the cavity as a function of t. A perfect coincidence then arises as two simultaneously 
detected photons at some time t. However, there are two challenges. One challenge 
is that there does not exist a perfectly efficient photodetector. Therefore, some pairs 
of photons are observed as single-photon emissions because one member of the pair 
escapes observation. In fact some pairs are missed altogether because both photons 
escape detection. The other challenge concerns the detection of two simultaneously 
created photons. Although created simultaneously, the emission from the cavity is not 
simultaneous due to the randomness of the emission time resulting from the nonzero 
cavity lifetime re -1 . 

We consider two photons to be coincident provided that they arrive within the 
window time interval r w . The choice of window time is not obvious, and it is our aim 
here to determine what the window time should be. As the two simultaneous photons 
can be separated by a time of order Ki . clS discussed above, the window time r w might 
be expected to be on the order of /c -1 . However, our purpose here is to consider the 
choice of r w in detail and to identify the optimal choice of window time r w that will 
produce the strongest 2PCR. 
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The choice of optimal window time is further complicated by the method of 
detecting nearly simultaneous photons. In the ideal case discussed above of a perfect 
photodetector yielding a record of all photon emissions from the cavity, one can then 
define a two-photon event as taking place if a second photon arrives between times to 
and t + r w , conditioned on a photodetection at time t . We refer to this rate as the 
conditional difference-2PCR and define this rate to be 

1 pto+T-w 

AS,(<J,^,r w )= lim - / dtA^(t,6,£ 1 ). (20) 

for A( 2 )(t,5,£i) = w^it^^Sx) - u/ 2 )(t,<5,£i = 0) with w®(t, 5, £ x ) given by Eq. (§. 
(This conditional difference-2PCR for a window time r w = was used in the quantum 
trajectory analysis of PCS in Carmichael et al. 1996). 

Another natural way to measure the 2PCR is by counting all photon pairs defined 
as being separated by an interval less than r w . This 2PCR is referred to as the 
unconditional difference-2PCR and does not rely on starting the count for the second 
photon conditioned on detecting the first photon. The definition of the unconditional 
2PCR is 

O pto+Tv, rt' 

A^ c (5,S u t w ) = lim — / dt' dtw®(t' (21) 

r 4 Jto Jto 

As shown in Appendix A , this expression can be simplified to read 

A^ c (6,£ 1 ,r w ) = - du (22) 

r w JO JO 

Both the conditional and unconditional difference-2PCR involve measuring photons 
and inferring whether two photons may be considered to have been 'simultaneously 
emitted'. This determination is based on whether two photons arrive within a specified 
time interval. In the conditional case, two photons are treated as being members of 
the same 'simultaneously emitted' pair if the second follows the first by less than the 
time window r w ; in the unconditional case, they are considered to be members of a 
pair if the temporal separation of the photodetections is less than the time window r w . 
Both definitions are reasonable and depend on the nature of the detection scheme. 
The window time for the conditional difference-2PCR will always be less than for the 
unconditional difference-2PCR. The conditional process is more efficient, in a sense, 
because the detection interval is always triggered by the first photon, whereas the latter 
case allows part of an interval to pass even before the first photon is detected. Thus, in 
the unconditional case, a longer interval is required to be able to detect the first photon 
and then subsequently wait for the second photodetection. 

We provide analytical solutions for the two extreme cases in |Appendix A| . The 
window time can be extremely long (kt w ^> 1), yielding expression ( |A.3| ), or extremely 
short (kt w <C 1), yielding expression ( |A.4j ) for both conditional and unconditional 2PCR. 
The short window time (r w — > 0) was the basis of the analysis of 2PCS in Sanders 
et al. 1997. In this treatment both the conditional and unconditional 2PCR at time t is 
approximated by (: n 2 (t) :). In the long-time limit the 2PCR is dominated by Poissonian 
statistics. 
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The choice of optimal window time r opt depends on the technique for observing 
two-photon coincidences, but another factor must also be considered. The purpose of 
2PCS is to observe two-photon decay resonances from the combined atom-cavity system. 
These 2PCR peaks are shown in Fig. |as a function of the normalized scanning field 
frequency 5. The choice of r w will depend on which peak is being observed. However, 
the best peak for observing a two-photon decay resonance occurs for 5 = 1 + V2. 

4.2. The Peak-to-Valley Ratio (PVR) 

In Fig. | the 2PCR peak at 5 = 1 + V2 resides over a background 2PCR that is largely 
independent of 5. In this Section we employ the parameters Si = £2 = k/2 (equality 
of Si and S2 helps to reduce the non-resonant two-photon absorption background as in 
Sanders et al. 1997) whereas, in Section |3|, Si = S 2 /2 = (as used in Carmichael 

et al. 1996). In this Section we choose lower driving field strengths which has the effect of 
producing smaller peak shifts than those studied in Sec. For this choice of parameters, 
the 2PCR background is (: n 2 :) ~ 2.1 x 10~ 5 . Let us characterise the quality of a 2PCR 
peak by the ratio of the peak height to the height of the background 2PCR. We can 
understand this figure of merit in terms of the ratio of signal to noise, where signal is the 
2PCR from the desired two-photon decay events, and the background noise corresponds 
to two-photon decays arising from unwanted off-resonance two-quantum excitations and 
decay events. The peak-to-valley ratio (PVR) is determined by the height of the peak 
to the height of the background (or valley) 2PCR. The optimal window time r w = r opt 
is defined such that the PVR for this 2PCR is maximal. That is, either a larger or a 
smaller choice of the window time would reduce the value of the PVR, thereby making 
the peak more difficult to discern from background events. 

There are other concerns besides the PVR in choosing the window time. For 
example, choosing a much shorter window time could improve the PVR but also lengthen 
the run time of the experiment required to accumulate enough signal. In other words, 
the absolute height of the peak is also a matter of concern in determining the feasibility 
of the experiment and is determined by the allowable timescale of the experiment. The 
minimum height would need to be on the order of T~ l for T the timescale of the data 
collection. 

The PVR is obtained numerically. The matrix continued-fraction method is used 
to solve the master equation to determine the peak height in Sanders et al. 1997. 
The background, or valley, is solved analytically by assuming a large detuning for the 
scanning field. The details are provided in appendix |Appendix B| . 

The 2PCR for large 5 is given by Expressions (|BT2[) and ( grp . The PVR for 2PCR 
is thus 

Af ] (5,Si,T w ) 

PVR 5 = * 2 w) (23) 
A^(r w ) 

with £ e {con, unc}. This convenient notation enables us to discuss both conditional 
and unconditional 2PCR with a single compact notation. In Fig. |9| surface plots of the 
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PVR vs g and r w reveal that the PVR increases as g decreases. This decrease is due 
to the background 2PCR (for 5 large) becoming small as shown in Fig. ^ whereas the 
pair signal rate does not decline as quickly for decreasing g. The increase in PVR is 
due to the fact that the resonant frequency for the transition |0) < — > |1)_ is u> — g, for 
a particular subensemble with coupling strength g, whereas the pump field frequency 
is constrained to satisfy uj\ = uj — gf. Hence, as g decreases, the pump field drives 
the system more and more off resonance. Thus, the likelihood of driving the system 
from |0) to followed by a photoemission, then repeating the driving |0) < — > |1)_ 
and a second photoemission, decreases rapidly with diminishing g: hence the background 
2PCR due to rapid driving to from |0) is quite small. A lesser contribution to the 
2PCR, due to 2ujx photons causing off-resonant excitation |0) « — > |2)_, also diminishes 
rapidly for decreasing g because the excitation moves progressively further from the 
resonance condition. Although the 2PCR due to de-excitation from |2)± decreases with 
diminishing g, the 2PCR for the background falls quicker thereby increasing the PVR. 

As our objective is to optimise the PVR, we choose the value of r w which, for each g, 
maximises the PVR. A single peak for each g is evident in Figs. |9](a) and (b). In Fig. [10 
we plot the value of the dimensionless quantity kt w vs g which maximises the PVR. 
This choice of r w would provide the optimal PVR for 2PCS if g could be fixed. It is 
clear from Fig. [10] that the window time r w for achieving the optimal PVR is an order 
of magnitude smaller than n" 1 . The rapid decline of r op t in the vicinity of g = gf = 9n 
is due to the driving from |0) to |1)_, followed by a photoemission, and repeating the 
process (as discussed above). Thus r w must be substantially reduced for g gf in order 
to excise the single-photon contributions to the measured 2PCR. 

Fig. [10] presents KT opt vs g/n for a range of values 0.2 < j/k < 10.0. The effect 
of varying j/k is to change the homogeneous linewidth particularly for the |1)_ level, 
which is important for the spurious contribution to the 2PCR, as discussed above. The 
rapid decline in r op t, as g approaches gf, varies with the choice of 7. For g small, the 
excitation from |0) to |1)_ by an u\ photon is sufficiently detuned that varying the 
homogeneous linewidth does not have much of the effect, hence the approximately flat 
curve for nr opt vs g/n in Fig. |lJJ However, the point of decline for «T opt does depend 
on the homogeneous linewidth, and we observe that the decline occurs for smaller g 
with j/k large, and the decline occurs for large g as j/k decreases. Fig. [l^ supports 
the earlier suggestion that the rapid oscillation between |0) and with concomitant 
photoemissions, is indeed responsible for diminishing the PVR. 

Fig. [10] is useful for understanding r opt , but the experimentally relevant choice of r opt 
is determined by taking the trace with respect to p in Eq. (|J). Using the same P(g) as 
throughout this paper, we obtain the PVR for the density matrix p. The highest PVR is 
plotted against Kr opt in Fig. [TT]. Fig. |TTJ depicts quite clearly the optimal choice r opt which 
maximises the PVR, and, for 7//C = 2, we observe that KT opt = 0.111 for conditional 
difference-2PCR and KT opt = 0.135 for unconditional difference-2PCR. 

Using this approach we plot KT opt vs 7/re in Fig. [TJ|. As P(g) is heavily weighted 
in favor of low g, the high-g decline in Fig. [TTJ is not so important. Hence, we observe 
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that r opt is relatively insensitive to the choice of 7, and the dependence is linear with 
correlation coefficients of —0.9983 and —0.9995 for the conditional and unconditional 
cases respectively. The dependence of r opt on 7 is given by KT opt = —1.4 x 10~ 3 7/k+0.11 
for the conditional difference-2PCR and KT opt = —2.1 x 10~ 3 7//t + 0.14 for the 
unconditional difference-2PCR. The variation of Kr opt for the range of 7/K is very 
small. It is clear that r opt is an order of magnitude smaller than Kf" 1 , and the 
dominant deleterious contribution to the difference-2PCR is the cycling |0) < — > |1)_ 
with corresponding photoemissions. 

5. Conclusions 

Photon coincidence spectroscopy (Carmichael et al. 1996, Sanders et al. 1997, 
Horvath et al. 1999, Horvath and Sanders 2001) is a promising technique for 
observing unambiguous experimental signatures of quantum field effects in quantum 
electrodynamics. This quantum field signature consists of two-photon spectral peaks at 
certain locations (with respect to detuning between the two chromatic components of 
the driving field). However, the peaks are not located exactly as predicted in previous 
studies due to multiphoton contributions to the two-photon count rate, due to non- 
simultaneity of the photon pairs emitted from the cavity and due to the effect of the 
cavity lifetime. 

We have considered these effects here. The multiphoton contributions are 
responsible for additional structure in the two-photon count rate peaks as shown in 
Fig. |8]. We observe that the higher-order contributions enhance the peak heights, rather 
than diminish them. Although additional structure is predicted by Fig. |^, for the chosen 
parameters, the additional structure is negligible on the two-photon count rate peaks. 
The technique of 2PCS is robust against multi-photon contributions to the two-photon 
count rate. 

Another important quantity is the detection window time r w . This quantity is 
the duration over which two photons are treated as though they originate from the 
same intracavity photon pair. For a separation of greater than r w , the photons are 
regarded as having been created independently. We consider two detection protocols, 
conditional and unconditional difference-2PCR. In each case the optimal window time 
is an order of magnitude smaller than the cavity lifetime due primarily to the need to 
reduce the |0) < — > |1)_ contribution to the two-photon count rate. We also note the 
negligible effect of varying the inhibited spontaneous emission rate 7 for the coupling 
strength distribution P(g) which is heavily weighted in favour of low g. 

These effects, namely multiphoton contributions to the pair count rate and the 
optimal choice of window time, are important in the design of photon coincidence 
spectroscopy schemes and in interpreting the results. Photon coincidence spectroscopy 
has been shown here to be quite robust, even with these potential deleterious effects 
being included in the simulations. 
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Figure 1. Two-photon excitation scheme from the ground state |0) to the first two 
excited couplets \n) (n £ {1, 2}, e G {— , +}) of the dressed states. The inhomogeneous 
broadening of the couplets |1) and |2) is 2% max and 2v / 2%max, respectively. Two 
two-photon excitations to the second couplet are depicted for a Dichromatic driving 
field with one component of amplitude E\ and the other with amplitude £2- The 
excitation pathway on the right (loi then L02) excites resonantly from |0) to |1)_ and 
then may excite resonantly to either |2)±. The excitation pathway on the left (u>2 
then uji) excites resonantly from |0) to |1)+ to |2)_ for g = — l)gj. 
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Figure 2. The 2PCR vs the normalized detuning S with £\/k = 0.5, £2/1 = 0.5, 
gf/n — 9, j/k = 2 for the inhomogeneously broadened system: w^ 2 \5,£i) is the solid 
line and A^(S : £i) is the dashed line. 
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Figure 3. Three-photon count rate (3PCR) vs normalized scanning frequency 
for g = gj = 9k, = 1/V2, SijrL — \/2 and j/k = 2. The solid line applies 

to (a^ 3 a 3 ), and the dotted line is for (a+a' 2 a 2 a-) . Alternating solid and dotted 
vertical lines are placed at values of S for which the 3PCR peaks are expected, 
namely, 5 G -(v^- 1), ~(VS- 1)/2, 2- y/Z, 1/V3, (V3+l)/2, V2+1, 2 + 
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Figure 4. The rate (a^ 3 a 3 ), which is depicted in Fig. ||, is reproduced here as a 
solid line in each of the four figures (a),(b),(c),(d). The rate (a^ 3 a 3 ) is repeated 
as a dotted line with £ x = for the |0) < — > transition in (a), with £2 = 

for |0) < — ► |1)_ in (b), with £\ — for the |0) < — ► |1)+ transition in (c) and 
with £ 2 = for the |0) < — ► |1)+ transition in (d). We present the four figures by 
allocating the two rows to each of the two transitions |0) < — > |1)± and the columns 
to whether £\ = or £ 2 — for the transition in question. 
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Figure 5. As for Fig. [|, with the solid line depicting the rate (a^ 3 a 3 ) and the dotted 
line corresponding to the same rate but with either £ x = (first column) or £ 2 = 
(second column) for the four transitions |1) E < — ► \2)' e with e, s' E { — ,+}. 
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Figure 7. Excitations to |3)_ for 8 = — 1, 5 = — l/v2 and for 6 = — l/v3 are 
presented here. For 5 = — 1, resonant excitation occurs to |1)_, followed by off- 
resonant excitation to |2)_, with a detuning of (2 — \2)gf = 0.59(7/, followed by 
off-resonant excitation to |3)_ with a detuning of (3 — \/3)<7/ = l-3p/. For 5 = — l/y2, 
off-resonant excitation from |0) to with a detuning of (1 — l/\/2)<7/ = 0.29g/, 

is followed by an resonant excitation to |2)_, followed by an off-resonant excitation 
to |3)_ with a detuning of (3/v2 — V3)gj = 0.39<?/. For 5 = — 1/V3, off-resonant 
excitation from |0) to with a detuning (1 — l/\/3)<?/ = 0.42(7/, is followed by an 
off-resonant excitation to |2)_ with a detuning of (v2 — 2/v3)fif = 0.26<7/, followed 
by a resonant excitation to |3)_. 
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Figure 8. Difference-2PCR vs normalized scanning frequency for an inhomogeneously 
broadened system and £\/n= l/v2 ; £2/1^= V2, j/k — 2. The solid line corresponds 
to A( 2 ) + A( 3 ) and the dotted line to A^ 2) . 





j/K KT ff X 10 J 

Figure 9. The pcak-to-valley ratio of the (a) conditional (PVR con ) and 
(b) unconditional (PVR unc ) 2PCR over the scaled coupling strength g/n and the scaled 
window time kt w for the scaled loss rate 7/fc = 2. 
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Figure 10. The scaled optimal window time KT opt vs the scaled coupling strength 
g/n for o j/k = 0.2, + j/k = 2, □ ~f/n = 5, x j/k = 7 and A j/k = 10 for (a) the 
conditional and (b) the unconditional 2PCR. 
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Figure 11. The peak-to- valley ratio (PVR) of the conditional (solid line) and 
unconditional (dashed line) 2PCR (for the masked atomic beam) vs the scaled window 
time kt w for the scaled loss rate r y/n = 2. 
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Figure 12. The scaled optimal window time KT opt vs the scaled loss rate j/k for the 
masked atomic beam with inhomogeneous broadening. The symbol + corresponds to 
the conditional r op t and o to the unconditional r op t- Linear regression methods yield 
the two lines. For the conditional case, the slope is —1.4 x 10~ 3 , the intercept is 0.11, 
and the correlation coefficient is r — —0.9983. For the unconditional case, the slope 
is —2.1 x 10~ 3 , the intercept is 0.14, and the correlation coefficient is r = —0.9995. 
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Appendix A. Conditional and unconditional two— photon count rate 

In the long-time limit, the conditional two-photon count rate (2PCR) is given by 

AgLO^i.T-w) = -(2kY / dt(: n(0)n(t) :> (A.1) 

with the trace being taken over p in Eq. (||). If the time window r w is large, compared 
to k -1 (the cavity lifetime), the two photons are highly decorrelated, and we can 
approximate 

(: n(0)n(f) :) — <n(0)> 2 . (A.2) 

Thus, 

Acon(^ £iT w ) — ► (2k) 2 (n(0)} 2 . (A.3) 

This count rate reflects the Poissonian nature of the count statistics for long window 
times. On the other hand, for kt w <C 1, the correlation between photon pairs cannot be 
neglected. Hence, the count rate reduces to 

A^ n (8,£ 1 ,r^^(2 K f(:n\0):), (A.4) 

which is the approximation employed in Sanders et al. 1997. 

Similarly, in the long-time limit, the unconditional 2PCR is 



2 



A^ c (S,£ 1 ,t w ) = -(2k) 2 alt' dt (: n(t)n(t') :) . (A.5) 



r 2 

This expression can be simplified as we show below. First we make the substitution u± = 
(t' ± t)/ v^2. We also introduce the notation d 2 u = du^du + and let V be the union of 
the two regions {0 < w_ < r w /v2, < u + < u_} and {r w /\^2 < u_ < v2t w , < u + < 
V2r w -u-} 

This substitution transforms the above double integral into the sum of two double 
integrals: 

AW( < y,£ 1 ,r w ) = 4(2«) 2 / / d\{:h{u'_)h{u' + ):) 
= 4( 2 -) 2 



d 2 u /: n(0)n(V2u + ) :\ (A.6) 



for u'± = (u+ ± tt_)/v2- The advantage of this expression is that the two-time photon 
number correlation depends on only one term in the double integral instead of both 
terms in the double integral. 

Greater simplification is possible and desirable for studying the short and long 
window time r w . Substituting u± = w±/\/2 transforms Eq. ( |A.6| ) to 



Aj»( tfj £ ljTw )=4(2«) a [ r*w- [ 

r w L^o Jo 



"2t w /■2r w — w— 

+ / dw- / dw 



dw+ 

(: h(0)n(w + ) :) , (A.7) 



o 
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which reduces to 

f*7 w ru 



9 /"Tw ru 

A^ c (g,5,r w ) = -(2K) 2 / du dw 

T w JO JO 



'0 JO 

x (: n(0)h(w) :) . (Ai 



For large (r w ^> n x ) and small (r w </t x ) window times A^nc reduces identically 
to Acoli as shown in equations ( |A.3|) and (|A.4j) . 



Appendix B. Background of conditional and unconditional 2PCR 

For the scanning field far off resonance (5 large), the time-dependent component of the 
Liouvillean £ can be ignored in the rotating picture. Thus, in the rotating picture, 
the master equation can be written as p = Cp with C, time-independent. The coupling 
strength g is fixed quantity, and averaging over P{g) (the inhomogeneous broadening 
case) is not considered in this appendix. If p is expressed as a vector, then £ 
can be expressed as a complex matrix with {— X n \n G Z N i\ the set of eigenvalues 
for N the number of levels in the Jaynes-Cummings ladder retained after truncation, 
and Re(A n ) > 0. The density matrix can be approximated by the sum 

N 2 

p(t) = Y,Pne- K{t - to) (B.l) 

n=l 

for {p n } a set of time-independent N x N matrices. Thus, the conditional 2PCR ([Z(]|) 
can be written as 

N 2 

1 f " ^ 

A ocon( r w) = c H / dt x^c n exp [-X n t] 

Tw J ° n=l 

N 2 

= c +^^{l-e-^}, (B.2) 

n=l 

with p n = X n T w and A n ^ for n > 0. Here the subscript on the left-hand side 
of Eq. ( p.2j ) is used to designate that we are considering the case of a monochromatic 
driving field. In the long window time limit, 

c = lim a£L(t w ). (B.3) 

T w >00 

Expansion (p.2j ) provides a useful method for calculating A{, 2 c ' on (r w ) by diagonalising 
the Liouvillean superoperator for the monochromatically-driven case. The function 
A 0con (r w ) is monotonically increasing because <9A 0con /cV w > if r w — > oo. Thus, 

/rj\ ____ 

dA Q l on /dr w — > as the function approaches the limit given by Eq. ( p.3| ). 
In the same way, the unconditional 2PCR can be obtained: 

A oLc( T ™) = Co + ^2 / du I dw 



^0 



N 2 

X 

n=l 



^ c n exp [-X n w] 
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In the r w — >■ oo limit, AQ„ nc (r w ) — >■ c . This is the same value for the background 
difference-2PCR as obtained for the conditional 2PCR. In summary, the Liouvillean 
superoperator is diagonalised. The master equation is solved in the t w — ► oo limit to 
obtain p ss = p(t — > oo). Then c = (2K) 2 Tr(a^ 2 a 2 p ss ). This quantity is the background 
2PCR for both the conditional and unconditional difference-2PCRs. 
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